The simplest quantum model supporting the Kibble-Zurek mechanism of topological 
defect production: Landau-Zener transitions from a new perspective 
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It is shown that dynamics of the Landau-Zener model can be accurately described in terms of 
the Kibble-Zurek theory of the topological defect production in nonequilibrium phase transitions. 
The simplest quantum model exhibiting the Kibble-Zurek mechanism is presented. A new intuitive 
description of Landau-Zener dynamics is found. 

PACS numbers: 03.65.-w,03.75.Lm,32.80.Bx,05.70.Fh 



5h 

xi 



T3 

S3 
O 

> 

O 
O 



^1- 

o 

-I— > 



i 

-a 
c 

o 



X 



In this Letter we present a successful combination of 
the Kibble-Zurek (KZ) [J Q theory of topological de- 
fect production and quantum theory of the Landau-Zener 
(LZ) model 0|. Both theories play a prominent role in 
contemporary physics. The KZ theory predicts produc- 
tion of topological defects (vortices, strings) in the course 
of nonequilibrium phase transitions. This prediction ap- 
plies to phase transitions in liquid He and 3 He, liquid 
crystals, superconductors, ultracold atoms in optical lat- 
tices 0, || , and even to cosmological phase transitions 
in the early Universe Q, Q|. The Landau-Zener theory 
has even broader applications. It has already become a 
standard tool in quantum optics, atomic and molecular 
physics, and solid state physics. The list of important 
physical systems governed by the LZ model grows. For 
instance, recent investigations point out that the small- 
est quantum magnets, Fe§ clusters cooled below 0.36K, 
are successfully described by the LZ model |(J. 

This Letter constructs the simplest quantum model 
whose dynamics remarkably resembles dynamics of topo- 
logical defect production in nonequilibrium second order 
phase transitions. The model is built on the basis of LZ 
theory and allows us to study the KZ mechanism of topo- 
logical defect production in a truly quantum case. Such 
a quantum insight into the KZ theory was up to now in- 
accessible except for the recent study of KZ theory in op- 
tical lattices filled with ultracold atoms In addition, 
we present a simple, intuitive, and accurate description 
of LZ model dynamics. 

For the rest of the Letter it is essential to introduce 
briefly the KZ theory. Consider a pressure quench that 
drives liquid 4 He from a normal phase to a superfluid one 
at a finite rate. Suppose the transition point is crossed 
at time t = 0, while time evolution starts at t <C 0. As 
long as the liquid is far away from the transition point 
its time evolution is adiabatic. In other words, the relax- 
ation time scale r, which tells how much time the sys- 
tem needs to adjust to new thermodynamic conditions, is 
small enough. As the transition is approached the critical 
slowing down occurs, i. e. r — > oo, so that at the instant 
—t the system leaves adiabatic regime and enters an im- 
pulse one where its state is effectively frozen- see Fig. ^ 
for an illustration of these concepts. The time t is called 
the freeze-out time and was introduced by Zurek 0. As 



the quench proceeds after crossing the transition point, 
the relaxation time scale decreases. At the instant t, the 
system goes back into an adiabatic regime. The freeze- 
out time is determined by the Zurek's equation: r(t) = t 
@. For the case of liquid 4 He it was found experimen- 
tally that t — To/|e|, where tq is a constant, while e is 
called the relative temperature. The latter measures the 
distance of the liquid from a transition point being at 
e = 0, i. e. e(t = 0) = 0. Physically changes of pres- 
sure translate into changes of e. It is further assumed 
that pressure changes are such that e = t/rQ, where tq 
is a quench timescale. Now the Zurek's equation reads: 
T oi~Q/i = t, which results in i — ^JtqTq. As shown in 0, 
knowledge of t allows for making a prediction of density 
of topological defects, resulting from a nonequilibrium 
phase transition, without solving dynamical equations 
describing the system! 

We consider time dependent Hamiltonian 
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written in the basis of time independent states |1) 
|2). Eigenstates of have the form 



(1) 



and 
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I (*)> 



cos(0(t)/2) sin(0(t)/2) 
-sin(0(t)/2) cos(0(t)/2) 
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where cos(0) = e/y/i + £ 2 , sin(0) = l/yl + £ 2 , 9 & 
[0,7r], e = A ■ t/uiQ. As in LZ theory A, w > are con- 
stant parameters. The level structure of Q is depicted in 
lower part of Fig. ^ while the gap equals y/ujfi + (A • t) 2 . 

Topological defects can be introduced into the LZ 
model in the following way. Suppose the state |1) corre- 
sponds to a vortex state being an eigenstate of angular 
momentum operator: £ z |l) = n\l) (n — ±1,±2, ...), 
while the state \2) satisfies L z \2) = 0. System's 
wave function can be written as |W) = a\l) + b\2) 
(|a| 2 + \b\ 2 = 1, — 8 l j). We propose to identify the 
density of topological defects with the average value of 



angular momentum: (\I>|.L z |\f r ) 



jjj. For the rest 



of discussion we define normalized to unity density of 
defects as Hi 
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Suppose now that the system undergoes adiabatic time 
evolution from the ground state of at t — ► — oo to the 
ground state of QJ at t — > oo. Therefore, the state of the 
system undergoes the "phase transition" from |1) to |2), 
i. e. from a vortex-defected "phase" to a vortex- free one. 
If time evolution fails to be adiabatic, which is usually 
the case, the final state of the system is a superposition of 
states |1) and \2) so that the final density of topological 
defects becomes non-zero. We will show that the KZ-likc 
theory predicts surprisingly correctly vortex density 
as a function of a transition rate only. 

(a) (b) 
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FIG. 1: Plot (a): relaxation time scale r from the KZ theory. 
Plot (b): inverse of the gap in the LZ model Q. Lower plot: 
energy levels of the Hamiltonian JQ; dotted line: cjq = case. 

Analogs of relaxation timescale, relative temperature 
and quench timescale are identified as follows. First of all, 
let us recall ingenious simplification of system's dynam- 
ics in the KZ theory. The simplification relies on the as- 
sumption that the system either evolves adiabatically, i. 
e. adjusts perfectly to changes of parameters, or becomes 
immobilized, i. c. undergoes the so called impulse evolu- 
tion 0- As proposed by Zurek, the switch between adi- 
abatic/impulse regimes is determined by the relaxation 
time scale, which is small during an adiabatic evolution 
and large in the impulse part. We would like to em- 
ploy similar strategy below. From the adiabatic theorem 
one knows that as long as the inverse of the gap is small 
enough the system starting evolution from a ground state 
remains in the ground state. It naturally suggests that 
inverse of the gap, being necessarily small in the adiabatic 
part of evolution, can be considered as a quantum me- 
chanical equivalent of the relaxation timescale introduced 
above: we set r = + (A • t) 2 |j. The equivalent 

of the relative temperature e, i. e. a dimensionless dis- 
tance of the system from anticrossing, is A • t/u>Q. As a 
quench timescale tq we take ujq/A, while loq we identify 
with 1/tq. Finally, we arrive at 

t=-== , £=—. 3 
VI + e 2 tq 

For |e| 3> 1 expressions (0 are identical as those intro- 
duced above in the context of topological defect produc- 
tion in liquid 4 He, which will be commented below. 



In the following we consider dynamics of the LZ model 
described by the Schrodinger equation: igjI'F) = H^S), 
with H given by We assume that time evolution 

starts from a ground state of Q at some t = ti and 
lasts till tf — * +oo. The quantity of interest will be den- 
sity of defects (J2J) at the end of time evolution, which is 
in fact the probability of finding the system in the ex- 
cited eigenstate at tf. Adopting the KZ simplification 
of the system's dynamics, we assume that the evolution 
of the system is either adiabatic or diabatic. The adi- 
abatic part takes place when the system is away from 
the anti-crossing, while the diabatic part takes place in 
the neighborhood of an anticrossing, where the inverse 
of the gap is so large that the system no longer adjusts 
to the changes of the Hamiltonian (Fig. ^>)- compare 
to the pressure induced quench in He described above. 
Therefore, the two non-trivial schemes can be considered: 

A: ti < —t: the evolution starts in the adiabatic 
regime, so it is adiabatic from ti till — t, then im- 
pulse from —t to i, and finally adiabatic from t to 
tf- see Fig. QJj, 

B: ti £ [—t,t]: the evolution starts in the impulse 
regime, therefore it is impulse from ti to t and then 
adiabatic from t to tf - compare to Fig. 

The statement that the evolution is impulse means that 
the system's wave function changes by the overall phase 
factor only. Additionally, it is assumed that we do not 
consider slow time evolutions for which the system stays 
whole time in the adiabatic regime due to the finite gap 
of the LZ model. Notice that such evolutions would be 
incompatible with KZ considerations where divergence of 
relaxation time scale at the transition point, Fig. [^i, pro- 
hibits adiabatic evolutions close to the transition point. 
The assumptions standing behind A (B) scheme classi- 
fication are approximate and heuristic as the whole KZ 
theory is, and our aim is to find how good they work in 
the LZ system. 

The only quantity that is still unknown is the instant 
i. It is found from the equation originally proposed by 
Zurek in the context of classical phase transitions 

T(i)=ai, (4) 

and modified by us by a factor a — 0(1), i. e. the only 
free parameter of our theory independent of tq and To- 
The solution of l@J reads 



e = e(i) 




The first observation shows that for fast transitions, i. 
e. x a —>■ at ro being fixed, one gets t = ^Jtqtq/oi. 
Therefore we recover, up to 0(1) factor, the well known 
result . It happens because in the fast transition limit 
i I and then r(tf) ~ tq/e, which is the same as in the 
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FIG. 2: Density of defects for the system whose evolution 
starts at the anticrossing center. Solid line- prediction 1101 . 
dots- numerical data. The parameter a = 0.77 was found 
from fit of il 1 ( )l) to numerics. 

theory of dynamics of quantum phase transitions in liquid 
4 He 0. This observation further supports similarities of 
our model to KZ systems. 

For the first application of our theory we consider 
the situation when time evolution starts far away from 
the anticrossing - a generic A scheme case. Taking 
= I I (*«)) as an initial system's wave function, 
and assuming limits ti — > — oo and tf — > oo, one easily 
gets the following final density of topological defects 

V n = \m f )\l}\ 2 « |<T (*)| I (-i))\ 2 = T^|2- ( 6 ) 
Derivation of © uses the following relations: 

IU(*/)I*(*/)>I « IU(*)I*(*)>I » KT (*)!*(-*)> I « 

KT (*)| i Substitution of I© into @ gives 



V » = ^T^ > V(x a )=xl+x a ^xl+A + 2. (7) 
V(x a ) 

Expanding T> n into a series one gets for fast transitions 
V n = exp(-x a ) + O (x 3 a ) , (8) 

which is an exact result up to (x^) terms if the con- 
stant a is chosen as it/2 Q- Notice that a = 0(1) as 
assumed in @. In the adiabatic limit (x a — > oo), Eq. 
predicts P n = 0(\/x\) instead of exponential de- 
cay, which does not affect results much due to very small 
value of D n in that regime. 

The best performance for fast transitions can be un- 
derstood as follows. The derivation of © requires as- 
sumption that in the time interval [— t, t ] the state of the 
system does not change essentially. The smaller is this 
time interval the better is this assumption. From © and 
|J5J) one easily finds that 1/tq grows monotonically with 
x a . Indeed, 1/tq equals ^fx^/a for x a — > and increases 
to 1/cv for x a —> oo- Therefore it is not surprising that 
our predictions work better for fast transitions. 

Now we would like to discuss the situation when time 
evolution starts from a ground state at the anticrossing 
center, ti — 0, which is a generic B scheme situation. As 
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FIG. 3: Density of defects. Solid lines- numerics, dashed lines- 
Eq. 11 U with a determined from fit. Upper (lower) curves 
correspond to 8 — 0.6%, a = 1.06 (do = 0.4-7T, a = 0.58). 



tf — * oo one gets 

V n = \{Ui)\im 2 = \(^-^=)i (9) 

where we put |*(0)) = | | (0)) = + ^|2), 

and assumed that |(T (*/)|*(*/)>| ~ KT (*)l*(*))l ~ 
|(t (t)|*(0))|. Combining © and © one gets 

V n = \ (l - v/1 - 2/V(x a )) , (10) 

with x a and V(x a ) defined in JSJ and 0. The agreement 
of this expression with results of numerical calculations is 
remarkable as depicted in Fig. [21 It is even better than 
in the previous case when we considered the evolution 
starting far away from the avoided crossing. We attribute 
it to the fact that now the frozen part takes less time, i. e. 
t instead of 2t, and to the absence of the approximation 
that the initial stage of evolution is adiabatic. 

We can also easily calculate density of defects when 
time evolution starts in the impulse regime (B scheme), 
but outside the avoided crossing center. Taking ^(ii)) = 
| i (U)) := -sin(0 o /2)|l) + cos(0 o /2)|2) we obtained 

cos(0q) 1-^1-2/7^(^)8^) 
v/27^) 2 ' 1 ) 

where 9o = arctan(cjo/(A - ti)) <E [0, 7r] measures distance 
of the starting point of time evolution from an avoided 
crossing, e. g. 9q — tt/2 when evolution starts from an 
anticrossing center and then Eq. (jll|l is the same as Eq. 

Comparison of (|llfl to numerics for tq/tq < 1.75 and 
\8q — tt/2\ < 7T/T0 reveals satisfactory agreement- see 
Fig. |3 for a typical situation. For larger tq/t and/or 
1 6*o — /2 1 the agreement gradually decreases, which we at- 
tribute to the fact that for these parameters the starting 
time moment, ti = tq/ tan0o, might be outside [— t, t], so 
that the assumption that the initial stage of time evolu- 
tion is impulse can be wrong. One avoids these problems 
when cither ti <C —t or \ti\ <C t, i. e. when the system 
evolves clearly within the A or B scheme, respectively. 

Having at hand above obtained results, let us comment 
on the the Zurek-like equality |Q extensively used in this 
paper. This equality gives the time moments ±1(tq), 
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which separate the adiabatic and impulse regimes |10| . 
Do we need to rely on this equation? To answer this 
question we notice that we aim at getting the best de- 
scription of LZ model dynamics by using the simplifi- 
cation that the evolution is either adiabatic or impulse 
in the sense specified below the A (B) scheme descrip- 
tion. It means that the whole problem can be reduced 
to getting the time moments ±t that lead to the best 
comparison of defect density to exact results. This can 
be done without Eq. by fitting the t(rg) directly to 
numerical data. Eq. reduces the problem of getting 
optimal t(ro) to fitting just one parameter, a, to the nu- 
merical results instead of fitting a whole function i(rg). 
Naturally, it is of fundamental interest to try to find ana- 
lytically the best t(jQ) function and compare it to those 
taken directly from the KZ theory. 

Let us point out possible usefulness of these consid- 
erations in understanding of non-equilibrium quantum 
phase transitions (QPTs). There is no doubt that quan- 
tum many-body systems of interest (spin systems, cold 
atoms in optical lattices) are more complicated than the 
LZ model. It is, however, generally accepted that close to 
the generic second order QPT point there exists an anti- 
crossing between a ground state and a first excited state- 
see Ch. 1.1 of P3 and Sec. 2.4 of 0. Therefore, one can 
expect that at least a qualitative picture of the change 
of system's properties during a QPT can be captured by 
a simple two level model, e. g. the LZ model. Notice 
that if it would happen that another sort of a two level 
approximation would work better, still it is quite likely 
that the same analysis as the one presented here would 
work. Once the many-body model of interest is specified 
and its low energy static properties are determined, one 
can define an equivalent of the "density of defects" and 
study system's dynamics using the tools presented in this 
Letter. In fact, the work, along the lines of this paper, 



on a dynamics of the superfluid-Mott insulator QPT of 
cold atoms in an optical lattice is in progress now. 

Let's look at other possible extensions of this work. 
First of all, it is desired to re-analyze more strictly the LZ 
dynamics, to get a systematic control of the intuitive re- 
sults obtained above. This work can be done utilizing the 
results from Second, it seems to be very interesting 
to investigate how the adiabatic/impulse simplification 
of system dynamics works in other quantum mechanical 
systems, for instance those which exhibit faster increase 
of the gap with the distance from anti-crossing. We ex- 
pect to get better qualitative agreement in these cases. 

Several other remarks are in order. First of all, we have 
shown that the very simple LZ model successfully repro- 
duces KZ-like dependence of topological defect density on 
the quench rate. Second, our results provide intuitive de- 
scription of LZ model dynamics unexplored up to now in 
numerous papers devoted to the LZ model. Notice that 
with our intuitive approach we have been able to obtain 
qualitatively correct predictions concerning LZ model dy- 
namics without solving the time-dependent Schrodinger 
equation! Third, we show that the KZ theory can pro- 
vide predictions beyond the lowest non-trivial tq/tq or- 
der usually considered 0,0, 01- Fourth, our results are 
directly applicable to different quantum two level sys- 
tems, e. g. to the molecular magnets Fes @- I n par- 
ticular, we expect that our prediction that the transi- 
tion probability for fast quenches (small tq/tq) decays 
in a power-like manner instead of an exponential one, 
Eq. (|10fl . can be experimentally verified and might be 
helpful in interpretation of experimental data. 
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